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ON  STABILITY  CRITBUA  OF  THE  KELLEY-MCSHANE  LINEARIZED 
THEORY  OF  YAWING  MOTION 


ABSTRACT 

A slightly  more  accurate  dynamic  stability  criterion  arising 
from  a more  careful  treatment  of  a pertuzbrtion  term  and  incorporating 
a revised  treatment  of  the  gravity  term  is  derived.  The  new 
criterion  has  the  property  of  reducing  td  the  classical  ballistic 
stability  requirement  s /when  only  is  retained  while  the  older 

criteria  became  indeterminate  for  this  case.  The  criterion  is  then 
gt;  eralized  by  the  insertion  of  the  requirement  of  a minimum  level 
of  damping  which  ie  not  neceeearily  zero.  Comparisons  are  made 
with  the  other  criteria  and  definition  of  the  terms  ••dynamic  stability1* 
"gyroscopic  stability"  and  "static  stability"  are  stated  and  dis- 
'*uesed.  The  effect  of  spin  on  dynamic  stability  is  indicated  and  the 
effect  of  the  location  of  the  center  of  mass  id  discussed  in  some 
detail.  A relation  for  an  optimum  center  of  mass  location  is  obtained. 
The  equations  of  yawing  motion  are  derived  in  an  appendix  and  solved 
in  the  body  of  the  report.  The  KelleyMcShane  force  system  is 
described  and  the  ^suiting  dependence  of  the  aerodynamic  coefficients 
on  center  of  mass  position  ie  derived  in  a second  appendix.  A third 
appendix  converts  the  expressions  used  in  the  body  of  the  report  from 
ballistic  to  aerodynamio  nomenclature. 
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INTRODUCTION 


The  Kelley-McShane  linearized  equations  of  yawing  .notion  are 
simple  linear  second  order  equations  which  possess  only  one  com- 
plication. This  is  caused  by  the  appearance  of  v,  the  spin  in 
radians  per  caliber,  in  certain  of  the  coefficients.  Since  this 
spin  is  a slowly  varying  function  of  position,  one  is  forced  to 
consider  the  effect  of  slowly  varying  coefficients.  In  order  to 
handle  this  Kelley  and  McSbane  J [6J1  introduced  a perturbation 
term.  The  effect  of  tills  tern  on  t»ie  stability  ha3  not  been  com- 
pletely discussed.  In  fact  most  of  it  is  omitted  from  tt ■ stability 
criterion  of  £57  and  some  it  is  omitted  from  that  of  £6]  2,  The 
clarification  of  this  discrepancy  and  the  resulting  deriva .ion  of 
a slightly  more  accurate  stability  criterion  is  one  of  the  aims  of 
this  report.  This  criterion  will  have  the  fu^+her  advantage  of 
incorporating  the  correct  contribution  of  gr*  .ty  as  is  indicated 
in  [2hJ 

It  is  found  that  west  considerations  of  stability  are  com- 
plicated by  unnecessarily  complex  subdivisions  of  stability  which 
arise  from  the  sign  of  the  moment  coefficient  and  the  presence  or 
absence  of  spin.  It  is  interesting  to  note  that  the  Kelley-McShane 
dynamic  stability  criterion  actually  becomes  indeterminate  when  the 
effect  of  all  aerodynamic  coefficients  other  than  K^  is  neglected. 

It  will  be  shown  in  the  body  of  this  report  that  the  Kelley-McShane 
criterion  can  be  revised  so  that  it  will  reduce  to  the  expected 
result  of  the  "classical”  ballistic  stability  requirement  s 
for  this  special  case  and  is  independent  of  any  restriction  on 
or  v* 


The  dynamic  stability  criterion  can  then  be  generalized  by  the 
introduction  of  the  concept  of  a minimum  level  of  stability.  The 
classical  stability  definition  io  founded  on  the  requirement  that 
the  two  damping  coefficients  and  be  non-n->gative.  The 

generalized  requirement  is  derived ‘from  the  requirement  that  and 

a2  be  not  less  than  an  assigned  number  a which  can  be  described  as 

the  minimum  permissible  level  of  damping. 


1. 

2. 


Numbers  in  brackets  refer  to  specific  publications  listed 
in  the  references  page£28J. 

To  be  precise  in  l?  I Sp  Bg  and  s^  which  appear  in  the  stability 


criterion  such  roceiva  a - Jg  from  the  perturbation  term  and 
in  [«J  they  receive  a - Jg  - md.  from  the  perturbation  term. 
J is  defined  to  be  the  " - gdu2  ^sin  0"  of  £5j  and  £6?. 
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A further  confusion  in  stability  studies  is  found  in  the  cany 
types  of  stability  which  are  discussed.  The  r.ost  common  of  these 
are  "dynamic  stability",  "static  stability"  and  "gyroscopic  stability*. 
Precise  definitions  of  these  terms  and  a discussion  of  their  in- 
terrelationship would  be  of  some  value. 


The  effects  of  spin  and  especially  of  center  of  mass  location  on 
dynamic  stability  of  a "statically  unstable"  missile  are  quite  in- 
teresting end  are_discussed  in  some  detail.  At  this  point  a dynamic 
stability  factor  s is  introduced. 

The  purposes  of  this  report  are  then; 

1.  To  derive  a dynamic  stability  criterion  which  contains  the 
full  contribution  of  the  perturbation  term  and  has  the  correct  con- 
tribution of  gravity. 

2.  To  state  this  stability  criterion  in  a fashion  which  Is 
independent  of  the  sign  of  1L,  and  the  presence  or  absence  of  non- 
sero  spin  and  which  is  equivalent  to  the  "classical"  ballistic 
stability  when  only  is  retained. 

3.  To  generalise  the  dynamic  stability  criterion  by  introduction 
of  a minimum  permissible  level  of  damping. 

U.  To  define  the  terms  "dynamic  stability",  "static  stability" 
and  "gyroscopic  stability*  with  respect  to  this  criterion  so  that  all 
interrelationships  are  dear. 

$•  To  disouss  the  effect  of  spin  and  center  of  mass  location  on 
dynamic  stability. 

SOLUTIONS  OF  THE  EQUATIONS  OF  YAWING  MOTION 


The  equations  of  yawing  motion,  which  ard  derived  in  the  appendix 
with  attention  to  proper  treatment  of  the  gravity  terms,  may  be  stated 

mat 

(1)  X"  ♦ ♦ Jg  - lvj  X»  ♦ - ivTj  X - G 

(2)  v*  - (b  - JT)  v 

where  X ■ Xj  ♦ iX^  complex  yaw 


md 


H - JL  - JD  ♦ £g-  ^ 


g 


i-  gd  sin  8 

V v 


1.'  fbi  s solution  is  essentially  that  described  in  the  original 
Kelley-Mc Shane  report  [$J  . 
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e 

V 


inclination  of  the  trajectory 


A 

•B 


M * md*  iL  * 

nr  15 

n 

I- JL- 


- - • nd 

"V  T 


M 


o-r*  - [ty>  - ->g  - T %>  * *]  r 

a * mass  of  the  missile 


d ■ diameter  of  the  missile 
B * transverse  moment  of  inertia 
A * axial  moment  of  inertia 

(g,,  g2,  g,)  » acceleration  due  to  gravity  vector  where 
x d * the  one  axis  is  along  missile's  axis 


“l 

*1 

r • 


• axial  component  of  velocity 

* axial  component  of  angular  velocity 

ir*2  ■ *i  v + ^*3  ■ *i  s3 d 

“ T“^  _ 


p ■ density  of  air 

are  the  ballistic  coefficients  and  primes  refer  to 

differentiation  with  respect  to  the  axial  arc  length  along  the 
trajectory  in  calibers. 


If  the  homogeneous  equation  is  first  considered,  the  problem 
reduces  to  the  solution  of  a second . order  linear  differential  equation 
* with  slowly  varying  coefficients.  Equation  (1)  may  be  simplified  by 

the  substitution  X ■ q exp 

the  first  derivative  term. 


1 /2j  -(H  ♦ Jg  - lv)  dp  which  eliminates 
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If  the  ? ic  considerations  are  employed,  equation  (1)  can  be 

written  {.0  >'  - x^q  - 0 

where  r2  - l/h  * 'rd&  - v2)  ♦ (H  ♦ J )2  *2  iv  (2T  - H - 

■V'  o ^ 


Equatio  * v-*/  baa  t>e  further  simplified  by  use  of  a second  transformation 

* 

w - (In  if)  . 

(U)  k*  * - - r2  ■ 0 


Clearly  *:  f r were  a constant,  a pair  of  solutions  of  (U)  is 
v *■  + r Gio*^/  examination  would  reveal  that  these  are  singular 

solutions  .-f  the  non-linear  equation  (It)*  The  solutions  of  the 
linear  (1)  which  correspond  to  these  are  linearly  independent  and 
hence  a linear  combination  of  them  would  be  the  general  solution  of 
that  equation.  Here  we  have  the  interesting  occurence  of  the  singular 
solution  of  one  equation  forming  the  general  solution  of  an  equation 
derived  from  ii..  Ve  now  assume  that  for  non-constant  r one  solution 
of  (U)  if  of;  »»hc  for*  r -3/2  £.  *&ere  £ is  small.! 

($)  r*  -1/2 E#-  £r  ♦ l/kf?  - 0 

r Q 

or  if  -1/2  ' * 1/U£*  is  small  in  comparison  with  r*  and  r^r  0 


<«  £-  f 

But  by  use  of  (2)  we  have 

(7)  -^7-  ' <D  * V 

r 2(lir2)  8 


fv  - i (2T  - H - D)  1 

S - 2iv  (2T  - H - DjJ 


i - v2  - hM  - (H  + Jg)2 


1*  The  -1/2  is  introduced  here  in  order  to  obtain  a slight  simplifi- 
cation in  the  final  form  of  the  solution*  This  perturbation  method 
is  called  the  W.K.B.  method  (Wentzel,*  Kramers,  Brillouin)  and  da 
described  by  K.  Jeffreys  in  Proc*  of  the  London  Hath.  Soc*  (2)  vol  23 
(1923)  p.  Ii28«  . 

2*  v -ljM  is  usually  written  as  v a where  a « 1 - - and  s,  the 

-2  s * 

ballistic  stability  factor,  is  • The  use  of  a squared  symbol 

implies  that  v2-  2|M  must  be  positive  and  the  definition  of  o2  has 
the  further  disadvantage  of  becoming  infinite  for  zero  spin*  In 
order  to  avoid  these  handicaps  m is  introduced* 
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V,  - 


H 


?-  second  solution  £f  the  fora  -r  -1/2C.  is  considered, 
tt  folio -6  <i Ircctly  that  £ • £ and  hence  two  solutions  of  (U)  are 

- -r  -l/2  if  . From  this  th&  g«nex*al  solution  of  (1) 


^ i,/  ft.  w ^3* 

is  obtaineu  as 


rp  f r r nV*J 

(8)  X • exp  1/2  j *j-H  - Jg  ♦ iv  - <-»  ♦ |_-n  ♦ 2iv  (2T  - H - D J J dp 

♦ *2  erp  1/2  J'33  - Jg  ♦ i7  - <5-  ♦ 2iv  (2T  - H - D^jdp 

If  we  write  £,in  the  fora  (L  *»  (£^  ♦ 1<F2),  ^ i3  c^ily  seen 

that  only^  should  be  consider  *d.  this  follows  froa  the  fact  that 
(D  - ,J_j  £ „i*  0*11  when  compared  with  v and  certainly  can  be  'I'-Usrui 
in  equation  (8). 

Now  (?)  £.,-  ffiLt j-L?I  , i - »2  -1^!  ^-(H  + J)2 

1 a2  ♦ It  v2  (2T  - H - D)2  ^ ^ ^ 

For  most  missiles  /» | r’|2T-H-0  j and  "v  < 1/20  and  we  have  the  usual 
T v2 

approximation  • 

X * 

Assuming  that  E,  T,  and  0 are  constants  it  can  be  shown  that  for 
a flat  trajectory  (Jg  - 0)  (8)  can  be  approximated  quite  accurately  by 

(ID)  X - ^ exp  (-c^  ♦ 10^)  p ♦ Kg  exp  (-«2  ♦ i02»)  p 
where  0^*  and  02*  linear  functions  of  p 
and  Og  are  constant  functions  of  p 

v - 0,*  09*  _ 


H 


<®1  ♦ ®2^-  “ .D 


M * 0-^’  • 0g * “ Og  ♦ D(o^  ♦ flg  — "Y" 

0i  - at 

-1/2  1(0,  - 0 ) -A 1 

D ■ . ttt— rr-s^r  (ses  note  1) 


[(«!- 

*1*  * >2 

*1 


• *2*  n 

: — 5—  -h-dJ 

« ♦ cL*  -/ 


• MV 


TT 


’ Itoe  &ret  four  relations  following  (10)  are  derived  fbr  the  case  of 
constant  0^*  * s'  in  £2!>J  • It  can  be  easily  shown  that  they  are  true 

for  the  case  of  linear  0^'  *s.  The  fifth  relation  conies  from  the 


.first  by  the  use  of  (2).  The  relation  fop  M 
j25/ due  to  the  retention  of  (h  * Jg)2  in  m.  5 


differs  from  that  in 


Although  the  solution  of  the  inhomogeneous  equation  is  not 
required  for  the  stability  discussion,  we  will  state  it  here  for 
completeness.  This  solution  adds  an  almost  constant  "yaw  of  repcse" 
Xpj  co  the  general  solution.  The  relation  defining  is  derived  in 

and  may  be  stated  in  the  form 

(u)  v J7jhs[: * ♦ 1 (w* ^ **Sf» 


-l 


where  0 is  the  inclination  of  the  trajectory. 

STABILm  CRITERION 


Dynamic  Stability  is  defined  by  the  statement  that  a missile  is 
dynamically  stable  if  the  yaw  described.. by  the  solution  o f the 

tnm  firtt.  H.f.nHn.  ts • 


homogeneous  equation  does  not  increase.  Referring  to  (8)  it  can  be 
seen  that  this  requires  that  H ♦ J ♦ (D  - JjiCi  ^ greater  than 

P D ^ ^ 


or  equal  to  the  real  part  of  the  square  root  and  be  non-negative* 
This  is  equivalent  to  the  relations 


(12) 


H ♦?«  j »£[-*  .♦  2iv  (2T 


- «- 


where  £•  J * (D  - Jg)£ ^ • Jg  (1  - £j)  ♦ • 

If  the  square  root  is  replaced  by  ^ a +ib  where  a - -in,  b - 2v  (2T  - H - D), 
a ♦ ib  "/(a*  ♦ b ; (cos  C ♦ i sin  C)  and  tan  C • <«  , 


De  Moivre 1 s Theorem  states  that: 

(13)  n({rm,)-*( 

■ui 


a2  ♦ b2 


(cos  | ♦ i sin  •£ 


i 


a2  ♦ b2 


c 

oos  -g 


Substituting  this  in  (12)  an$  using  the  half  angle  formula  for  cos 
gives  the  following  equationst 


(1U) 

H 

or 

h 

or 

H *£$■ 

a2  ♦ b2 


J 


1 + cos  C 

£ 


* b 


1. 


2a 


Since  the  initial  conditions  affect  only  the  homogeneous  part  of 
the  solution,  this  requirement  la  equivalent  to  the  restriction  that 
the  size  of  yaw  caused  by  the  initial  conditions  will  not  increase 
throughout  the  trajectory.  Although  the  yaw  of  repose  may  increase, 
it'  is  unaffected  by  the  random  initial  conditions  and  its  effect 
may  be  computed. 

Note  that  £ is  here  replaced  by  (D  - 0£v 

D 





*3*2222232221 
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cola 


Squaring  xhi.a  we  hare 

(15)  2 fH  -t-Sj  2 • -j£2  + b2  ♦ a,  H ♦ £ 2 0 

The  second  inequality  is  required  in  order  that  (15)  be  equivalent  to 
(Hi).  Regrouping  and  squaring  again  we  h awe 

(16)  [2(H  ♦C)2  -a]  2 * a2  ♦ b2,  H ♦£-  0 

Unfortunately  (16)  is  not  precisely  equivalent  to  (15)  since  the 
first  inequality  of  (16)  implies  either 

(17)  2(H  +t  )2  -a  * ♦ b2  or  a -2(H  *S)2  *7 fiT*  b2 

In  order,  that  (16)  be  equivalent  to  (15) , it  must  be  possible  to 
rejeot  tho  second  inequality  of  (l?).  But  this  Is  oily  possible 
if  H +tfO.  If  H ♦£  * 0#  (15)  requires  that  b • 0 and  that  a 
be  non-positive*  With  this  in  mind  we  can  revise  (16)  to  the 
following  fora  which  is  equivalent  to  (15)  • 

(16»)  H ♦!>()'  H 0 

b - 0 

or-  <0 
J5(h  ♦!)*-. 2 .b* 

If  the  inequalities  are  solved  for  -a,  and  a and  b are  replaced  by 
their  definitions, 

(Ifl)  H «£>0  H ♦IT m 0 

or  v(2T-H-D)  - 0 

- e2  -I|M  " 0 

v2  >2iM  ^ fkQsJL~?ll2-  (h  *£)2  ♦(H+J  )2 

L H ♦ £ J * 

These  are  the  exact  stability  conditions.  The  usual  size  approximations 

ban  now  be  applied  which  state  J2  terns  can  be  dropped1  and  the  follow- 
ing theorem  is  obtained* 


1*  Iota  that  for  constant  spin  D - * 0 or  aero  spin  v ■ 0, 

these  J2  terns  vanish  exactly  from  (16).  Hence  the  theorem 
is  not  an  approximation  of  (18)  for  this  case* 
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Theorem  In  the  linearized  theory  a missile  posses ses  dynamic 
stability  if  either  of  the  following  sets  of  relations  is  satisfied 


(l?a)  H +Fr  0 


\ h 


(19b)  H +£  • 0 

v(2T-H-D)  • 0 


vhsre  £ * Jg  (l-£^)  4 

£ - vl  *C(w2  -1<H)  ♦ 2(2T-H-D)2] 
1 (v2  — ItM)2  ♦ Uv^(2T-H-D)2 


v -W9  0 


-2 

v 


V2  -liM  . 


v2 -IjM  > 1 2T-H-D  j 

Note.  (19a)  nay  be  obtained  somewhat  more  simply  by  the  assumption  that 
t n|  « 

n » * -liM^O  together  with  the  approximation  of  a binomial  expansion 
of  the  square  root  in  (8).  If  this  is  done,  the  a ^ and  of  (30) 

become 


(20) 


0^2  - 3/2  [h  *Z  J =|r  (2T-H-dJ'<'  1/2  [(1  1 1)  H7  |t  ♦ 


(<{i2  • 1/2  » (X  1 O) 


rm~ 

where  <r  * <| w 

v 


p Dv2  . D 
m 


7" 


Dynamic  stability  requires  that  2 0 and  this  implies  that 

a * s i laagau 

. fT 

(21)  J vZ  - IjM  g [v(2T-H-D)l  « But  this  is  equivalent  to  (19a) 

H ♦ £ 

Although  this  derivation  is  handicapped  by  the  assumptions  made  above, 
with  some  modifications  this  reasoning  can  be  made  rigorous* 
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r mms 


f 


j 

i 


\ 

f 


* 


—2 

Since  relations  (19)  require  v -kN  to  be  positive,  £ ^ must 

be  positive  and  considered  as  a function  of  spin  it  attains  a 
minimum  of  one  for  infinite  spin.  If  the  slope,  D-J  , is 
1 — G 

positive  then  the  minimum  value  oft  ^ corresponds  to  the  minimum 

t and  the  following  sufficient  stability  condition  follows  from  (19a) . 


(22)  H ♦ D*  0 


o * v2f2T-H-D 

» ■**  Lh.dJ 


if » * jd  - js. 


is  approximated  in  “*»>»'  JD>  there 
results  a set  of  conditions  which  are  sufficient  for  stability  over 
flat  trajectories  when  • 

(23)  H ♦ JD>0 


v2 


Equations  (23)  are  precisely  the  stability  criteria  of  (£] 

with  the  correct  treatment  of  gravity?  and  equations  (22)  ara  uhe 
criteria  of  \6j»  Therefore  if  only  flat  trajectories  are  con- 
sidered and d ^ <.  , the  stability  criteria  of  [$J  and  £6j.are 

sufficient  conditions  for  stability?'  No  indication  of  hie  second 
set  of  relations  of  (19)  can  be  found  in  (Sr’*)  or  (?!)  and  herein 
is  the  reason  ibr  the  failure  of  the  criteria  of  nr  falto 
reduce  to  the  classical  gyroscopic  stability  crito.ua,  s £ 1 • 


1*  This  is  satisfied  over  all  of  the  upward  branch  and  part  of  the 
downward  branch  of  any  trajectory  and  is  true  over  the  entire 
length  of  flat  trajectories  dealt  with  in  spark  range  work*  The 
minimum  value  of  £ . is  one  and  hence  the  minimum  c is  D. 

2*  In  part  of  £ isTjategrated  out  and  no£  considered  in  the  stab- 
Itv  ur.  \lysis.  From  the  remainder  the  correct  treatment  of  gravity 
£2^) produces  (23)  which  is  independent  of  gravity. 

3*  It  must  be  emphasized  that  from  a practical  engineering  standpoint 
(19) > (22)  and  (23)  are  equivalent  and,  hence,  (22)  being  the 
simplest,  should  be  used* 
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before  going  on  to  a discussion  of  various  types  of  stability 
and  the  effect  of  spin  and  center  of  mass  location  on  stability 
we  will  first  obtain  a generalized  form  of  (19)  which  should  be 
quite  useful  in  design  work.  (19)  was  obtained  from  the  require- 
ment that  2^  °*  ***  wer®  replaced  by  the  requirement  that 

a.  a whefe  a is  the  minimum  permissible  damping,  the  criterion 
beiomes  much  more  flexible  and  valuable  to  the  design  ballistician 
or  aerodynamicist. 

This  revision  can  be  done  quite  easily  by  replacing  (12)  by 

(12')  can  be  reduced  to  (12)  by  writing  £ for • £ -a.  This 
means  that  the  generalized  criterion  can  be  obtained  by  inserting 
C-«  fbr  £ in  (20), 

Theorem  In  the  linearized  theory  the  damping  rate 3 of  a 
missile'  are  greater  than  a when  either  of  the  following  sets  of 

relationlia  satisfied.1 

(2l*a)  H ♦ l -a  > 0 

■ 

v2-l*$  ' 

L H +C-  a 

DISCUSSION  OF  STABILITY 

In  addition  to  the  concept  of  dynamic  stability  two  other  types 
of  stability  are  usually  considered:  gyroscopic'  and  static.  Their 

definitions  follow  directly  from  the  definition  df  dynamic  stability 
for  the  oase  where  only  the  effect  of  is  considered.  For  this 

casq  the  second  set  of  relations  of  (19)  must  be  used  and  they  reduce 
to  . , 

(25)  v2  -UM$  0 


(2l*b)  H ♦ L -o  - 0 
v(2T-H-D)  - 0 

v2  -1*M£0 


(12*)  H *C 


-H 


+ 2iv  (2T-H-D) 


1/2 


9 


1.  Note  that  a can  be  negative  although  its  magnitude  is  restricted 
to  the  order  of  a J term.  This  means  that  should  a little 
instability  be  permitted  the  generalized  criterion  will  still 
apply. 


H* 


M m U'WWjWHM.; 


A missile  is  said  to  be  gyroscopic-illy  stable  if  it  satisfies 
inequality  (25).  For  zero  spin  (25)  becomes  H or  JC,  4 0.  A 

missile  for  which  ds  0 is  said  to  be  statically  statue.  For 


statically  unstable  missiles  K.  ^ 0 and  inequality  (25)  can  be 

written  in  the  formt 


(25') 


V1 2 

TH 


2 2 

.A2 

UB  pd^  4 ^ 


This  is,  of  course,  the  classical  ballistic  stability  condition 


s 5?  1 where  a » — • Gyroscopic  stability  for  statically 

UB  pdJ  Uj*  K^j 

unstable  missiles,  therefore,  reduces  to  the  classical  stability 
condition. 

The  logical  connections  between  gyroscopic  stability,  static 
stability^  and  the  basically  fundamental  dynamic  stability  can  be 
derived  from  their  respective  definitions*  They  are  shown  in 
figure  1 and  are  listed  below. 

• 

(1)  If  only  Kw  is  considered  for  stability,  gyroscopic  stability 

is  sufficient  for  dynamic  stability.  ^ 

(2)  Gyroscopic  stability  is  always  necessary  for  dynamic  stability. 

(3)  If -spin  is  aero,  static  stability  is  necessary  for  gyroscopic 
stability. 

(U)  For  any  spin,  static  stability  is  sufficient  for  gyroscopic 
stability. 

(5)  If  H * 0 and  v ■ 0,  static  stability  is  sufficient 

far  dynamic  stability. 

(6)  If  H and  2T  - H - D ■ 0,  gyroscopic  stability  is 

2 

sufficient  for  dynamic  stability. 


1.  Since  Kjj  is  defined  so  that  it  is  positive  when  the  center  of 

pressure  is  in  front  of  the  center  of  mas3,  most  bodies  of  re- 
volution  are  statically  unstable  while  most  finned  missiles 
are  statically  stable. 

2.  This  relationship  is  quite  interesting  since  it  indicates  a sort 
of  optimum  center  of  mass  location  i.  e»  such  that  2T-A-D  • 0. 
This  relationship  is  implied  by  R.  Turetsky  in  [l8j. 
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DYNAMIC  STABILITY* 


(1)  With  damping 

H *Z  * 0 


(2) 


or 


IH  «■  J 5 0| 
g 


v ■ 0 


or 


GYROSCOPIC  STABILITY 


v2  -l|H  s 0 


v •»  0 


I 


JTATIC  STf'TUTr 


Jig,  1 


* S - jg  (1-  £3)  ♦D61 


£.  v2fe!  -UP  ♦ 2{ZML-1ifJ 

n (v2-  UM)2  * UvZ(2T-H«D)Z 


v2 


vfc  — Um 
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j.7.  WT  >rr*  w 


Without  damping 

h ♦ r - o 

v(2T-H-D)  - 0 
v2-!*  2 0 


F 

[ s 


In  order  to  discuss  the  particular  effect  cf  srin.  we  will  use  the 
slightly  simpler  sufficient  condition  (22)  and  will  assume  H ♦ D to 
bo  always  positive.  (22)  can  bo  written 

(26)  v2  ^ ^ M or  v2  a(2-s)^  i*M 

(H  ♦ D r 


- 2T 

where  s - 3-7-5 


2(K^  - k2  Ky) 

h * 3 % - 


(Dymrdc  Stability  Factor) 


k^  • -j-  (k^  is  axial  radius  of  gyration  in  calibers) 


*2* 


md 


■g-  (kg  is  transverse  radius  of  gyration  in  calibers) 

n 

If  s is  either  0 or  2*  the  coefficient  of  vc  is  zero  and  for  this 

case  a statically  unstable  missile  can  not  be  made  dynamically  stable 

by  spin  while  a statically  stable  missi.le  is  dynamically  stable 
regardless  of  spin.  ..Otherwise  (26)  can  be  solved  for  v and  two 
possibilities  ariset 


(27a)  a(2  - s)7  Oj 
• (27b)  s(2  - s)<  Oj 
From  this  we  see  that! 


v2* 


v2* 


UM 


"5(2  - "s) 

liM 

s(2  - s) 


(1)  For  a statically  stable  missile  either  "5(2_-  s)«$.  0 and  it 

is  dynamically  stable  independent  of  spin  or  s(2  - s)*  0 and  dynamic 

stability  places  an  upper  bound  on  spin  by  use  of  (27b).1 

(2)  For  a statically  unstable  missile  either  s(2  - s)>0  and 

dynamio  stability  places  a lower  bound  fbr  spin  by  use  of  (27a)  or 

s (2-  "a)  $ 0 and  the  missile  is  dynamically  unstable  independent  of 

spin.2 

1*  Dynamic  data  obtained 'from  the  simple  finned  configurations  de- 
scribed in  07}_show  that  this  upper  bound  actually  exists*  At  Mach 
number  of  1.7  s is  -7.2  and  for  this  value  of  i the  upper  bound  on 
spin  is  U»2°  per  caliber.  One  round  has  been  fired  at  a spin  rate 
of  U.0°  per  caliber  and  it  was  found  that  the  slow  arm  had  practi- 
cally no  damping.  At  Mach  No.  of  1.3,  s ■ - 30  and  model  was  unstable. 

2.  This  entire  discussion  of  the  effect  of  spin  on  stability  is 
essentially  that  of  Kelley-McShane  in  either  [5J  or  [6].  The 
S£  *s  of  Kelley-McShane  may  be  related  to  the  symbola  of  this 
report  by  s.  • H ♦ £ , Sg  " 2T  +£  -D,  and  " 2(H-T)  +€  * D. 

The  discussion  is  inserted  here  in  order  to  give  this  report 
completeness*  Here  again  the  case  of  the  airship  must  be  excluded. 

(See  Append*-:  B). 
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For  the  assumption  made  above  (ramely  H ♦ £ “ H + D 7 0], 
the  effect  of  spin  on  dynamic  stability  may  be  stated  in  a very 

concise  manner*  Theorem  If  s(I  - s)  and  M have  the  same  sign, 

* -2 

v*  > ■* 

a missile  is  dynamically  stable  if  s - — % 
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Otherwise  spin  has  no  effect  and  statically  stable  missiles  are  dy- 
namically stable  while  statically  uns  table  missiles  are  dynamically 
unstable* 


We  will  now  consider  the  interesting  question  of  the  effect  of 
center  of  mass  location  on  the  stability  of  a statically  unstable 
missile.  In  order  to  do  this  we  need  to  convert  the  generalised 
stability  condition  into  a more  useful  and  specialized  form. 
Straightforward  algebra  and  the  approximation  £ i D can  be  employed 
to  obtain  the  following  result. 

Theorem  , 

The  damping  exponents  ou  and  Ck  of  the  epicyclic  yawing  motion 
of  a statically  unstable  missile  are  greater  than  or  equal  to  an  , 
assigned  value  a if  either  of  the  following  sets  of  relations  are 
^ satisfied.  (If  0*0,  the  relation  of  the  first  set  became 
precisely  those  for  the  dynamic  stability  given  in  equation  (27*)! 


(28a)  H ♦ D - a > 0 


(28b) 


-2 

v 


wr 

sv2-a) 
(28c)  04  j<2 


M / . 2T  • a 

where  s - s (a)  - +■'  g"-'o 


or  (29a)  H + D - a - 0 

(29b)  2T  - a - 0 

(29c)  s&.l 

(Generalised  Dynamic  Stability  Factor) 


Relatione  (29)  Are  only  of  ACAdemic  interest  for  designers  And  will 
be  disregarded  in  the  remainder  of  this  report.  The  fora  of  relations 
(28)  has  been  greatly  simplified  by  the  introduction  of  ’5(a)  which 
will  be  called  "the  generalised  dynamic  stability  factor".  Whenever 
Is  appears  without  a reference  to  a value  of  a,  a will  be  taken  to 
be  sero  and  1 is  then  called  "the  dynamic  stability  factor".  If 
we  multiply  the  numerator  and  denominator  of  s(e)  by  -i--~-.it  reduces 


to  the  following  simple  form* 


ptb 


*0(0*)' 


2(1^  - ny2  up  -o 

■4.  ♦ *»■*.%  -*i'z  h - °* 


where  a"  ■ 


'2 

m 


Pd; 
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The  equation  a • — ■ is  plotted  in  Fig-arc  2.  The  curve 
s(2  - s) 

is  a fairly  flat  one  with  asynpbdtea  at  zero  and  two.  The  modification 
of  the  classical  gyroscopic  stability  requirement  by  the  dynamic 
stability  factor  is  significant  only  when  a is  outside  Jjhe  interval 
from  3A  to  !>A*  • The  center  of  this  interval  which  is  *5  ■ 1 
corresponds  to  the  optimum  center  of  mass  position  mentioned  before 
(2MI-D  ■ 0)  and  8 is  independent  of  a here.  This  means  that  if 
the  gyroscopic  stability  factor  is  unity  for  this  center  of  mass 
position,  an  increase  in  spin  which  in  turn  increases  ”s"  will  cot 
improve  the  dynamic  stability.  As  a natter  of  fact  equation  (20) 
shows  that  both  arms  damp 'at  the  saps  rate  which  is  independent 
of  spin  and  that  this  rate  isH  ♦ t H ♦ D.  Equ  ation  (20)  also 
reveals  the  interesting  property  that  for  141  the  slow  arm  has 

the  smaller  damping  rate  while  for  1 jrl  the  reverse  is  true.'1’ 

This  is  important  since  a of  the  criterion  is  the  damping  rate  of 

the  am  possessing  the  smaller  damping  rate. 

, ‘ . 1 

Va  now  return  to  the  center  of  mass  effect  on  stability.  From 
appendix  A we  have  the  result  that  and  are  independent  of 

and  Kg  are  linear  and  quadratic 

functions  respectively  of  its  location. 

(30)  T> 1Q  ♦ T1  q 

n ♦ D - Hq  ♦ I^q  ♦ H2  q2 
where  Tq  - JT 

T1  "^1  2 JF 

. ho  - ♦ y 8 ^ V ja 

H1  * V*  <J5  * 


1.  By  alow  or  fast  arm  is  meant  the  arm  .with  the  smaller  or  larger 
turning  rate  The  alow  am  is  sometimes  called  the  precesslonl 

arm  while  the  fast  arm  is  called  the  nutationd  arm.  Since  in  the 
classical  theory  of  the  top  the  processional  motion  is  usually  con- 
sidered to  be  the  motion  of  the  plane  of  yaw  and  the  nutational 
motion  is  the  variation  of  the  magnitude  of  yaw  fbom  maximum  to  minimum, 
this  nomenclature  is  conflicting.  (It  can  be  shown  that  nutation  rats 
is  actually  A*,-  (fa  while  the  processional  rate  is  variable  unless 
K.  or  K9  is  *>ro)* 


center  of  mass  location  while 


DYNAMIC  STABILITY  FACTOR 


q is  the  shift  of  c.m.  in  calibers  and  is  bounded  to  a variation  over 
which  the  radii  of  gyration  may  be  considered  constant. 

Relation  (27a)  now  defines  the  limits  on  q in  order  that  dynamic 
, stability  be  possible:  The  requirement  that  s(2-s)  be  positive  may 

be  written  as: 

♦ T-.  q , 

(31a)  0<s<2  or  0*  ± *— < I 

Ho  * H1  * + H2  q 

For  supersonic  velocities  and  most  bodies  of  revoluation,  ex- 
perience shows  that  for  *5  • 0 the  center  of  mass  is  to  the  rear 
of  the  centroid.  Beyond  this  point  which  is  not  more  than  two 
calibers  away  from  the  centroid  s < 0 and  it  is  impossible  to 
stabilize  by  spin.  The  situation  can  of  course  be  improved  by- 
changing  k, . For  some  long_bodies  of  revolution  (fineness 
ratio  greater  than  eight)  s • 2 when  the  center  01  mass  is  forward 
of  the  centroid*  Forward  of  this  point  s > 2 and  it  is  impossible 
to  stabilize  by  spin*  This  situation  can  be  improved  by  changing 
kg  an<}/or  1^* 

In  conclusion  we  will  construct  two  representative  stability  plots 
• for  long  bodies  of  revolution.  The  second  plot  ia  quite  similar  to 
those  obtained  for  shorter  models*  Table  I giver  values  of  the 
aerodynamic  coefficients  and  physical  constants  for  a typical  body 
• -and  Figure  3 plots  the  square  of  the  required  spin  against  c.m* 

location  for  various  values  of  a*,  v*  is  used  In  prder  that  s ■ 1 
will  appear  as  a straight  line.  A value  o*  spin  diich  is  above  an 
[ , a*  curve  corresponds  to  a missile  with  damping  exponents  greater  * 

i than  a. 

I ‘ ' TABLE  I 

1 . . 


*A  - >01 

Kp-.lS 

V -3 

«H-a> 

*M-1 

%•  -9 

it,.  - -.3 

CO 

-a 

CM 

V 

CM 

• 

■ 

CM 

1 

jf 

4 -S*io-S: 

Notice  that  there  is  only  a small  range  of  c.m*  for  which  dynamic 
stability  is  possible*  Point  0 is  the  optimum  c.m.  location  and  is 
about  *U  cal  rear  of  the  centroid*  If  we  now  attempt  to  Improve  the 
stability  characteristics  of  this  hypothetical  model  by  increasing 
•2 

kg  from  *2  to  *3,  we  find  that  the  situation  .is  completely  changed 

(see  Fig  U.).  The  optimum  point  is  .1  cal  forward  of  the  centroid 
and  the  region  of  possible  dynamic  stability  is  tremendously  enlarged. 
Note,  however,  that  the  restriction  on  c.m.  location  reappears  when 
damping  of  the  order  of  tt*  • 6 is  required. 
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The  major  effects  that  the  consideration  of  the  dynamic  .-♦•ability 
of  statically  unstable  missiles  has  on  the  development  of  such  spin 
stabilized  missiles  may  be  summarized  in  the  following  three  statements: 

(1)  If  there  is  no  requirement  on  the  size  of  damping  and  3/U^-s  < 5/U, 
the  usual  gyroscopic  considerations  apply  (s  Z 1). 

(2)  If  there  is  no  requirement  on  the  size  of  damping  and  "5  is 
outside  this  interval,  it  must  be  improved  either  by  a different  mass 
distribution  which  varies  k,  and  or  by  a different  shape  which 
changes  the  ballistic  coefficients. 

(3)  IT  there  is  a minimum  level  of  damping  a,  then  1s(a)  must  be 
considered  and  the  design  procedure  is  that  of  (l)  if  3/U  *•»(«)  <•  S/h 
otherwise  it  is  that  of  (2). 

C.tfr  Tru^yti^ 

C*  H.  MURPHY 
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TABLE  OF  SYMBOLS 


A Axial  moment  of  inertia 
B Transverse  moment  of  inertia 
C • arc  tan 

D'  Jp-¥-JA 
F ■ aerodynamic  force 

o-r  -£»„ - Jg - t *>  * <7 * 

H • (hphgjh^)  angular  momentum  vector  or 


Jt  ~ Jts  * 


md 


mu  * 

5“  n 


H0’H  jiHg  coefficients  in  expansion  of  H * D as  function  of  center  of  mass 


« “r 

Axial  spin  deceleration  coefficient 
Xjj  Drag  coefficient 

Axial  drag  coefficient 
Kp  Magnus1  force  due  to  yaw  coefficient 
Kjj  Damping  moment  coefficient 
Lift  coefficient 


1*  Terms  in  the  cross  force  (A2)  or  cross  moment  (AU)  which  are 
aero  for  aero  spin  are  called  "Magnus*  forces. 


2$ 
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K^,  Overturning  (righting)  moment  coefficient 

Kjj  Normal  force  coefficient 

Kg  Damping  force  coefficient 

K Magnus  moment  coefficient 
T 

/ 

K_  Magnus  force  due  to  cross  angular  velocity  coefficient 

Kjj,  Magnus  moment  due  to  cross  angular  velocity  coefficient 

Ki  1C,  Integration  constants  in  (10) 
x»  * 2 

M ' " HT  JM 

(M^jMgjM^)  Aerodynamic  moment  vector 
T * JL  • TT  JT 

T0*T  ^ Coefficients  in  scansion  of  T ae  function  of  center  of  mass 

a • - m 


b • 2»  (2T-H-D) 

d diameter  of  missile 
g acceleration  due  to  gravity 
radii  of  gyration 

(81*62**3)  sector  acceleration  due  to  gravity 
^\*h2'h3^  fcnsular  momentum  vector 
m mass  of  missile  " 

. r4  u, 

p * J dt  independent  variable 


q ■ X exp  l/2^*(H+Jg  - iv)dp  or  shift  in  center  of  mans  location 
r*  1/2  £(UM  - v 2)  ♦ 2 iv  (2T  - H - d£) 
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8 stability  factor 

✓ 

a dynamic  stability  factor 

SpSgjS^  symbol*  used  in  stability  discussion  of  [Sj  and  £6j 
t time 

u (u^UgjUj)  vector  velocity  of  center  of  mass 
v - (In  q)f 

a minimum  level  of  damping 
Opt-.  damping  of  each  epicycle  arm 

£ £ • W.K.B.  perturbation  term 

Cm 

6 inclination  of  trajectory  to  horizontal  plane 
Ug  u- 

X * — ♦ i complex  yaw 

T.  T. 

Xjj  Taw  of  repose 
eagd  « d 

|i  ■ rr—  ♦ i complex  angular  velocity 

T.  . T. 
cad 

v * — — spin  in  radians, per  caliber 

- T. 

v-  |v 

p density  of  air 

c symbol  used  in  f5J  and  frSj  see  footnote  to  Equation  (7) 

#\t  ftli  to  , to  turning  rates  and  their  derivatives  for  the 

x c A * epicycle  arms 

-O.  angular  velocity  of  coordinate  system 

(eopCOpCi^)  angular  velocity  of  missile 
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APPENDIX  A KELLEX-McSHANB  FORCE  SYSTEM 


In  order  to  describe  the  aerodynamic  force  acting  on  a missile 
in  flight,  we  need  first  to  set  up  a suitable  coordinate  system. 
This  will  consist  of  numbered  axes  so  that  the  1 axis  lies  along 
the  missile's  axis  while  the  2 and  3 axes  are  perpendicular  to  it. 

As  a convenient  shorthand  we  will  consider  the  plsne  perpendicular 
to  the  missile's  axis  to  be  the  complex  plane  with  the  2 axis  as 
the  real  axis  and  the  3 axis  imaginary,  now  if  the  Telocity  vector 
of  the  missile's  center  of  mass  is  resolved  along  these  exes,  the 
yaw,  X,  of  the  missile,  which  is  the  directed  angle  between  the 
missile's  ax^s  an^the  velocity  vector,  can  be  expressed  in  complex 

form  as  X • ^ — 1 ^ ■ . This  expression  is  of  course  accurate  for 

small  yaw  only. 


The  second  and  third  components  of  the  angular  velocity  vector 
<0  - (ooj,  09g,  ssj)  are  of  importance  in  the  consideration  of  the  yawing 

motion.  These  can  be  non-dimensionalised  and  written  in  complex  form 

09-  d ♦ i 09-  d . 

as  i*  * where  d is  the  diameter  of  the  projectile. 

The  aerodynamic  problem  can  now  be  reduced  to  a dynamic  problem  by 
the  assumption  that  the  aerodynamic  force  and  moment  are  functions 
of  the  air  density  p,  it*  sound  velocity  a,  axial  velocity  u. , 
axial  angular  velocity  oa^,  yaw  X,  cross  angular  velocity  p,  size 

of  missile  which  is  characterised  by  its  diameter  d and  its  shape. 

The  KelleyMc  Shane  theory  makes  the  further  simplification  that  the 
aerodynamic  ibree  and  moment  are  linear  functions  of  X and  ji.  If 
these  hypotheses  are  applied  to  a configuration  possessing  an  angle 
of  rotational  symmetry  lads  than  180°  and  a plane  of  mirror  symmetry, 
the  aerodynamic  force  and  moment  can  be  written  in  the  following  form 

m.  c«j 

(A1J  1A  * -pd2u£ 

(A2)  P2  4 1P3  • p d2  u2  Q-*  4 iv  Ky)  X 4 (v  4 iKg)  4] 
(A3)M1--pd3^vKA 

(AU)  Mg  4 1M3>  p d3  £(-v  Xj  - iXjj^X  ♦ 4 ivK^ 


where  v 


» 
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An  important  aspect  of  the  above  definitions  is  the  dependence 
of  the  coefficients  appearing  in  (A 2)  and^Ak)  on  the  location  of  the 
center  of  mass.  In  (A2)  the  aerodynamic  force  is  independent  of  the 
location  of  the  center  of  mass  while  the  yaw  X is  defined  in  terms 
of  the  velocity  of  the  center  of  mass.  The  dependence  of  the  moment 
coefficients  in  (AU)  is  more  complex  since  they  relate  to  the  total 
aerodynamic  moment  about  the  center  of  mass  and  are  also  associated 
with  X. 

Consider  the  case  where  the  center  of  mass  is  shifted  from  the 
point  0 on  the  missile* s axis  to  the  point  0*  which  is  located  on 
the  missile *8  axis  at  a distance  of  q calibers  from  0.  (Positive 
q will  correspond  to  a shift  toward  the  nose  of  the  missile.)  The 

vector  00*  is  in  coordinate  form  (q  d.  0,  0).  All  quantities  relating 
to  the  configuration  possessing  the  center  of  mass  will  be 
marked  by  an  asterisk. 

Now  if  corresponding  points  of  the  two  configurations  possess 
the  same  motion*  the  total  aerodynamic  force  on  each  configuration 
are  identical  and  the  total  aerodynamic  moment  when  computed  about 
corresponding  geometric  points  will  be  equal  [3} , [6j,  [26f  For 
the  same  motion  we  now  compute  the  velocitiesSi  ti  e ooint:;  0 and  0*. 

(A$)  u - u ♦ cox  (qd,  0,  0)  - (x^,  «•  qdj  u^  - <»2  qd) 

Therefore 

(A6)  ^ 

(A7)  X*  - X - iqn 

Since  the  angular  velocity  vector  ie  independent  of  the  location 
of  the  center  of  mass*  we  have 

(A8)  v*  » v 
(A9)  n*  • p. 

We  can  now  express  the  force  acting  on  the  configuration  with  the 
new  center  of  mass  as 

(A10)  - -p  d2  Uj*8  i -p  d2  Vj*  S^* 

(AU)  F*  ♦ lFj*  • pd2  u*  [(-*/  . iv*  K/)  A*  * (,*  Kjp*  ♦ Uj.*) 

- p d2  g-l^*  ♦ i»Kj* )(A  - lop)  ♦ (»  Kxf*  ♦ iKg*)  f) 
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If  we  equate  the  above  with  (Al)  and  (A2)  and  conpare  coefficients 
of  A and  p respectively  we  obtain  the  following: 


(Al 2) 

V 

■*» 

(A13) 

h* 

-*■  . 

(A1U) 

h* 

■h 

, • 

• 

(A15) 

h* 

-Kg-qX,, 

• 

(Al£) 

" *XF'  qKF 

$ 

Now  in 

order  to  compare 

the  aerodynamic  moments  we  must  first 

compute  the  moment  for  the  second  configuration  about  the  point  O. 

This  can  be  written  as  the  moment  about  0 in  terms  of  the  starred 
quantities: 

(A17)  M*  ♦ (<W/  0,0)X  F%  M*  - (qd,  0,  0)XF  - (h£,  M*  - F.J  qd,  - F 3 qd) 
Therefore 

(A18)  - -p  Uj2  d3  K* 

(A19)  - F3  qd)  ♦ i (M^  ♦ Fg  qd)  - M*  ♦ iM*  ♦ iqd  (Fg  ♦ iF3) 

- pd3  u^  ^(-vkj  - iK^)  ( X - iqn)  ♦ (-^  ♦ ivKxr!)  * 

- q [(-%  *'1»  Xp)  X ♦ (v  Kjj.  ♦ 1K-) 

Equating  (Alfl)  and  (A19)  with  (A3) and  (All)  we  have 


(A20) 

_* 

ka* 

** 

(A21) 

V* 

*M" 

*«- 

(A22) 

V* 

kt* 

itj  - 

qXy 

(A23) 

_* 

*H  * 

« {Ks  * V + **  *» 

(A2U) 

K * 
XT 

' XXT 

” 0 (*»  * *rJ  * •>*  *f 
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In  order  to  obtain  the  form  of  (Al)  and  (A2)  used  in  Appendix  B 
we  derive  two  Xinal  relations.  Sometimes  the  force  system  is  reduced 
to  only  two  components  which  lie  in  the  plane  of  yaw  and  are  resolved 
along  the  trajectory  and  perpendicular  to  it.  The  first  component 
is  called  the  drag  force  0 and  the  second  is  called  the  lift  force  L. 
They  define  two  mors  coefficients  by  the  relations 

(MS)  D - -p  d2  u,2 
(A26)  L - -p  d*  u£  ^ X 

v 

If  we  neglect  Magnus  terms  And  terms  dependent  on  fi  in  (A2),  ve 
see  that  (Al)  and  (A2)  describe  the  same  type  of  force  as  (A25)  and 
(A26)  but  differently  directed.  We  now  select  tie  real  axle  so  that 
it  lies  in  the  plane  of  yaw.  Clearly  (A25),  (A25)  .are  related  to 
(Al)  and  (A2)  by  a rotation  through  the  angle  X. 

• • (A2*»)  D ■ Fg  sin  X ♦ F^  cos  X 

(A26)  L - Fg  cos  X - F^  sin  X 
Hence  for  small  X,- 


. ...  (A27)  Kp-XjjXsinX  + Kj^cosXaK^ 

(M8)  \ - o°.  X - Vk,,  - Kj, 

If  the  discussion  of  the  above  paragraph  scorns  to  lack  rigor, 
equation  (A27)  and  (A28)  may  be  taken  as  definitions  of  Kq  and  X^. 

- This  can  be  done  since  only  these  relations  and  not  their  physical 
definitions  are  used  in  this  report.  Using  these  ve  can  write 
(Al)  and  (A2)  in  their  final  form: 

(Al*)  Fx  - -P  d2  u2  Kp1-. 

(A2*)  ?2  * iF^  ■ p d2  - Kp  ♦ v Kp)X  ♦ iXg)  t* 


33 


APPENDIX  B:  DERIVATION  OF  THE  EQUATION  OF  TAWING  MOTION 


A right  handed  orthogonal  coordinate  system  with  axes  mustered 
1,  2,  3 moving  with  .the  missile  and  so  orientated  that  the  1 axis 
always  points ’toward  the  nose  along  the  aissile1  s axes  will  be  used 
throughout  this  appendix.  If  we  specify  » * («j,  Og*  co^)  to  be  the 

angular  velocity  vector ^ of  the  aissile  in  this  coordinate  system 
and  take  JX"  U2i>J22'AP  **  th®  *nSnlar  velocity  of  the 

coordinate  ays  tea,  the  above  restriction  on  the  1 axis  provides  the 
- relations/)  2 • andJl^  - ay  is  selected  to  be  identically 

s arc  and  th*  2 axis  is  initially  orientated  to  lie  in  the  horizontal 
' plane  and  point  to  the  right*  For  the  small  yaw  and  small  curvature 
of  the  trajectory  which  are  assumed  in  this  development,  the  2 axis 
will  always  remain  quite  close  to  the  horizontal  plane*  Our  equations 
of  motion ; aTet  ’ ;; . ' . ‘/V. . 

(M)  (m  u)  - (Fj,  Fa,  Fj)  ♦ ■ (g^,  gj,  g3) 

(88)  % ’ 


where 


a ■ (u^,  Ug,  u^)  is  the  velocity  vector  of  the  center  of  gravity 

1 " *^3)  the  angular  momentum  vector  of  the  missile 

(F  , F0,  FO  is  the  aerodynamic  force  vector  acting  on  the 
1 * J missile 

(g^,  gg,  g^)  is,  the  acceleration  due  to  gravity  vector 

(Mi,  Na|  M3)  is  the  aerodynamic  moment  vector  acting  about  the 
' * ' ' center  of  gravity 

m is  the  mass  .of  the  missile  ... 

Sines  our  coordinate  system  is  not  an  inertia  system,  we  have  to 


1*  All  linear  and  angular  velocities  are  with  respect  to  an  inertia 
system  fixed  on  the  earth*'  • 


3U 


fmmmm 


”***■* *vur 
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differentiate  using  the  relation  d m d * Six. 

cFt  at  fixed  axes 

Although  axis  1 is  fixed  in  the  missile  the  other  two  are  moving  with 
respect  to  it.  This  would  make  the  angular  momentum  vector  rather 
complicated  were  it  not  for  the  fortunate  fact  that  the  mass  distribution 
of  most  missiles  are  effectively  those  of  bodies  of  revolution.  This 
implies  that  every  direction  normal  to  the  axis  of  the  missile  is  a 
principal  axis  of  inertia  and  one  transverse  moment  of  inertia  prevails. 
If  A is  the  axial  moment  of  inertia  and  B the  transverse  moment  of 
inertia,  then  (hp  hg,  h^)  is  the  vector  (A Mp  Bcog,  Boo^). 


It  is  convenient  to  make  the  plane  normal  to  the  missile’s  axis 
a complex  plane  with  axis  2 the  real  axis  and  axis  3 the  pure 
imaginary  axis.  For  any  vector  equation  this  can  be  done  by  multi- 
plying the  third  coordinate  equation  by  i and  adding  to  its  second. 

In  the  remaining  equations  of  this  section  a dot  will  signify* 
differentiation  with  respect  to  time* 


With  the  above  remarks  in  mind  it  is  possible  to  obtain  the 
following  equations  from  (Bl)  and  (B2). 


# F 

(B3)  * (e^  Ug)  - £l  + g 

. . , F9  ♦ i F, 

(BU)  (Ug  ♦iuj)  - i^  (« g ♦ i «^)  * ■ * (gg  ♦ i g^) 

* (B$) 


(B6)  (ig  ♦ i Wj)  - i («jg  -f  i ec^/  ^ 


Mg  ♦ i Mj 

— s 


2.  If  we  take  Xp  Xg,  Xj  to  be  a triad  of  unit  vectors  pointing 

along  the  respective  axes  of  the  coordinate  system,  any  vector 

V can  be  written  as  V ■ ♦ VgXg  •**  T^X^.  Differentiating 

with  respect  to  time  using  the .dot  convention  this  yields 

V ■ (fcjX^  ♦ VgTg  ♦ V^Xj)  * (Vi^l  ♦ VgXg  ♦ V^).  The  first  term 

and  to®  second  can  be  written  as  SL  X V.  This  is 

vSt/fixed  axes 

the  definition  for  the  angular  velocity  vector  o£  the  coordinate 
system.  In  component  fornJUs  (Xg  . X^)  X1  ♦ (l^  . X^Xg  ♦ . Xg)  » 


3$ 


The  equations  are  now  made  non-dimensional  by  use  of  the  axial 
velocity  u^  and  the  missile’s  diameter  d.  This  introduces  the  rcn- 

disnensional  dependent  variables* 


o>jd 

“ 


(spin  in  radians  per  caliber) 


tu  ♦ i 

X ■ .■  — ■■ — - (complex  yaw). 


COg  d ♦ i co^ 


(complex  angulAr  velocity) 


dt 


which 


The  independent  variable  tins  is  replaced  by  p - F *1 

o “T 

for  small  yaw  is  approximately  arc  length  measured  along  the  ° tra- 
Jeotoiy  in  calibers.  If  a prime  is  introduced  to  represent  differen- 
tiation with  respect  to  p,  we  have  the  relation  ( )*•  ( ) * ^ . From 

equations  (B3)  - (B6)  we  can  now  write* 

<w)i:.  .*id  -*id 

~~  "T- 

/ u'l  .(^2  4 * *3>  d (e2  4 i *0  d 

- (B8)  X*  ♦ —i  X - i*»  « —2 1 ♦ -J: 

*1  »«i2  V 

V VL  dZ 

(B9)  v*  * .4.  v - — i— ■— 

*1  A tu2 

H A (iu  ♦ i Mj  dZ 

l*  • 

Rrom  equations  (Al*)i  (A21),  (A3)  and  (Ai*)  the  linearised  force  system 
can  be  written  as*  . .....  . 

,1  > *«  V 


(BIO)  n« 


(bu)  - 

(B12)  F2  ♦ i F3  - 
(B13)  — — m v 


JD  ♦ iv  JF)X  ♦ (vJjjj,  ♦ i 


(Bill)  M2  + i M3  - n tu,2  Q-  v JT  - i JM)  X'+  (-JH  + iv  Jxr^j 

where  J •<  K for  all  subscripts,  p is  the  air  density,  and  the 

K's  are  the  aerodynamic  coefficients.  (See  Table  of  Symbols  and 
Coefficients) 


Using  (Bll)  we  can  now  write  for  (B7) 

(BJS)  - - JD  ♦ Jg  - Su3~“3ug)d 


where  J 


*ld 

m mAmm 

6 S7 


T 

»1 


For  the  small  yaw  and  transverse  angular  velocity  required  by 
the  linearity  assumption  the  third  term  on  the  right’  side  of  (61$) 
can  be  neglected  in  comparison  with  ~Jq'+ 

• ui 

. . (B16)  3--  . JDv  J 

a * . 

Using  (Bl6)  and.  definitions  (B12)  - (Bill)  equations  (B8)  - (BIO) 
may  be  rewritten  asi1  •. 

(B17)  X*  - i n • ( - ^ ♦ i v JF)  X ♦ (v  Jjp  ♦ i J3)|i  ♦ T 

(B18)  v»  - (D-Jg)  w 

(B19)  |»>  - (JB  - Jg  ♦ »)  |1  - IcJ2  jj-  v JT  - i JM)X  ♦ (-JJ,  ♦ lv  JH)nJ 

vtor.  y - <*2  * 1 «3)d  . X - " ®1X2)  + i (g3  - g1  X^J  d 

y—  2 


^ jA 
jj-2  r a / 


(k^  is  axial  radius  of  gyration  in  calibers) 

■ — g-—  (kg  is  transverse  radius  of  gyration  in  calibers) 


1«  The  grouping  of  J X in  the  y was  a correction  Introduced  by  Professor 

o 

McShane  in  order  to  make  the  treatment  of  gravity  more  accurate. 

A more  extended  discussion  of  this  may  be  found  in  £4 
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We  now  operate  on  (Bl?)  with  the  operator}  4r  ♦ ( j - J_  + J ) 

I up  c ti  u g 

• i » (|  ♦ kg2  multiply  (Bl?)  by  v J^p  ♦ i (1  ♦ Jg),  and  add. 

Assuming  that  derivatives  of  force  coefficients  can  be  neglected  and 
using  (B18)  to  eliminate  v*,  the  result  reduces  tot 

(B20)  X"  ♦ Jjl  ♦ Jg  - iv  xj  X*  ♦ (-M  - iv  T)  X - 0 ♦ v 

where 

H ■ h - jd  * .S*  h 

»■!» 

X • 1 ♦ l^2  JjJ  ♦ | Jy 

M * kj2^  * y2  kj2  Jp)  ♦ v2  UT  Jjj  - JT  Jff) 

+ J3  * -iO  * (2D  * V 

- tmh  -h*  •>*  4 '*1*  [*  * - JXP  ^ * JS  Jt  * h.  ,XX  * I VaJ 
o • r'  - [fy  - nj2  •$,  - rfg)  ♦ tv  (1  ♦ kj2  JjjH  t 


The  upper  case  letters  with  the  exception  of  Q are  selected  in 
order  to  identify  the  moment  coefficient  which  is  the  principal  con* 
stituent • the  quite  formidable  expressions  above  can  be 
by  certain  quite  reasonable  site  assumptions.  Ve  assume 


wwir 

i simplified  i 
that  j Jp  j <10"a, 

|jtJO  x ID*3  otherwise,  1<,^<102  , k^2<»*vO,  kjj24.  2,  | v|^l,  and  that 

|X  |<  1.  Since  £JL  i*  usually  about  $ x 10 this  restricts  the 

magnitude  of  Kp  to  less  than  2 and  that  of  other  K^'s  to  less  than  60» 

The  requirement  for  the  special  case  of  reduces  to  (10*2u) 

where  d is  in  feet  and  u in  feet/sec.  From  (B17)  it  can  be  seen  that 
p la  comparable  with  iX'  and  hence  the  second  term  on  right  side  of 
(B20)  can  be  neglected  in  comparison  with  iv  XX'.  X can  clearly  be* 
very  well  approximated  by  1 and  similar  approximations  apply  to  the 
other  terms.  There  results  the  following  good  approximation  of  (B20). 

(B21)  X*  ♦ (H 


whore 
H * 


JL*JD 


♦ J - i v)X'  ♦(-H-vT)X*Q 

, 1 < . 


'JM  * * 
-2  , 


T " " ^1 

0 - r’  - 0JD  -■  *5*  -to  • V * 1 'Jr 
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Note  that  the  fora  of  (B21)  differs  from  (B20)  only  by  the 
absence  of  X and  Since  the  general  stability  analysis 

of  the  report  is  stated  in  terms  of  M,  H,  and  T,  all  the  assumptions 
of  (B21)  except  X » 1 and  vJjy(D-Jg)ik  ■ 0 may  be  avoided  by  inserting 

the  values  of  H,  M and  T.  given  after  equation  (620). 


An  ijqxjrtant  case  where  the  above  assumptions  do  not  apply  is 

that  of  the  airship.  Here  p d^  is  of  the  order  of  one.  If  the 

* 

effect  of  drag  and  gravity  aVe  neglected  and  spin  is  taken  to  be 
aero  equation  (B20)  reduces  to: 

t 

<B22)  X"  ♦ HX»  - MV  - 0 
shore  H ■ 


and  Jg  is 
negative 

even  when  Kg  is  positive.  This  provides  the  interesting  result  that 

a statically  unstable  configuration  can  be  dynamically  stable 

without  spin^.  This  is,  of  course,  limited  to  configurations  of 
snail  density. 


Since  for  this  configuration  and  are  positive 
usually  negative  and  leas  than  one.  we  see  that  M can  be 


1.  An  explicit  example  of  this  is  given  in  Durand,  Aerodynamic  Theory. 
Vol.  VI.  pages  UO  - 112.  The  negleotion  of  drag  and  gravity  is 
quite  valid  since  drag  is  usually  neutralised  by  the  ship's 
. propulsion  system  and  gravity  by  the  ship's  buoyancy.  The 
assumptions  that  the  center  of  buoyancy  ia  located  at  the  center 
of  mass  and  that  apparent  mass  of  the  air  due  to  flow  around  the 
ship  oan  be  neglected  are,  however^  implicitly  implied  in  (B22). 
Since  the  center  of  buoyancy  is  usually  located  above  the  center 
of  mass  and  ths  additional  consideration  of  apparent  mass  affects 
M by  not  more  than  10  JC,  $22)  is  still  reasonably  valid  for  the 
horisontal  component  of  yaw,  1^.  Actually  only  the  behavior  of 

this  component  is  considered  in  Durand. 
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APPENDIX  Cs 


CONVERSION  FROM  BALLISTIC  TO  AERODYNAMIC 
NOMENCLATURE 


The  work  of  this  report  has  been  done  in  terras  of  the  ballistic 
K'e  which  are  little  known  outside  the  field  of  ballistics  and  may  be 
quite  confusing  to  an  aercdynamiciat  who  does  his  dynamic  stability 
analyses  in  terns  of  the  aerodynamic  C' s.  It  is  therefore  quite 
worthwhile  to  express  the  results  of  this  report  in  terms  of  these 
symbols.  This  effort  is  handicapped  however,  by  the  three  facts* 

1.  The  missiles  usually  treated  in  ballistics  have  a rotation 
symmetry  which  results  in  pairs  of  aerodynamic  coefficients 
being  equal  and  hence  corresponding  to  only  one  ballistic 
coefficient. 

* 

2.  In  ballistics  the  missiles  usually  have  a high  rate  of  spin 
and  Magnus  effects  have  to  be  considered  to  which  there  are 
no  corresponding  aerodynamic  coefficients. 

3.  Terms  involving  the  rate  of  change  of  angle  of  attack  appear 

in  most  aerodynamic  stability  analysis  while  no  such  terms 
appear  in  the  usual,  ballistic  force  system.  t , 

The  axial  conqxsnenta  of  the  aerodynamic  force  and  moment  are 
usually  defined  in  aerodynamic  nomenclature  asi 


«a) 


X-- 1/2.2*  SC, 

V-  1/2  P v2  S b($)  C, 


(02) 


From  this  we  see  that 

Kj,  - 1/2  S/d2  Cjj 

*A  • - VU  (sb2/dVc^ 

(pb/2V) 

If  the  transverse  components  of  the  aerodynamic  force  and  moment 
are  assumed  to  be  linear  functions  of  yaw,  change  in  yaw,  and  angular 
velocity,  and  Magnus  coupling  is  introduced,  we  have  the  following 
definitions* 


X - (1/2  p V2  ^ 

*[%. 


* 4 CXr  4 CY^ 


4CI 


pq 


W ♦ o^]  (®j 


hO 


(C3) 


(1/2)  p V2  S^£cz  a * CZ  ^ * CZ. 

*['y  ’ • Vs  ’ V*?®) 


M - (1/2  p V2  Sc)  |C  a ♦ C (■§£)  ♦ C.  'C° 


mn  mn 

a q 


“a 


*[\»  P * Sr  ^ 

H - (V2  P V2  S b)fjcap  P ♦ (|?) 

+ [cv  “ + * v^0 

If  the  missile  is  asswnd  to  porsrss  trigonal  or  greater  rotational 
symmetry,  it  follows  that  [h],  [6],  [*>],  [27 Jt 


(CU) 


\ ■ \ 

% ■ % “ v 

* % * CN, 

%,  '%J'  \q 

' % * \ 

v ■ %p  $ * V 

%♦  ' *\ 

%p  • v|}  * V 

* \ 

-c  (~)  - C (-)  a -CL 

V C V*  ^pq 

C (-)  ■ C (-)  jfc  CL 

la 


m 


The  third  aet  of  symbols  is  introduced  in  order  to  emphasize  the 
existence  of  symmetry  and  will  be  employed  throughout  the  remainder  of 
this  appendix.  If  we  insert  these  symbols  into  (C3),  multiply  the 
second  and  fourth  equations  by  i and  add  to  the  fii-o  and  third  respectively 
there  results: 

r ♦ iz  - (1/2  p v2  s )frcv  ♦ 1 cN  1(0  * **) 

u pq  q J 

HV‘« 

♦ 1* - d/z  p l2  esmfr)  0^  - 1 C^(p  ♦ to) 

If  equation  (CJ>)  is  conqpared  with  equations  (A2)  and  (Al*),  the  Magnus 
and  non-Magnus  static  coefficients  are  easily  related. 

...  ; ' •'«'  ...  , V ...••■■  ' ' , . . 

*U  • - 1/2  s/d2  Cj, 


« 1/2  Sc/d?  Cjj 


(C6)  - Xp  - 1A  Sb/d®  cN 

P® 

«*  ’ f »'  ‘V‘ .j'  1 ■ *1  V ■ ‘ 

■ ; • - lA  Scb/d1* 

The  relationships  between  the  remaining  dynamic  coefficients 
are  somewhat  more  complicated.  Fortunately  it  can  easily  be  shown 
that  the  remaining  Magnus  coefficients  are  lost  in  the  differential 
equations  of  yawing  motion  due  to/ths  J2  connection.  It  therefore, 
remains  only  to  connect  two  ballistic  coefficients,  Kg  and  K^with 


U2 


four  aerodynamic  coefficients,  Cjj  , 

q,  a q d 


In  order  to  do  this  we  need  only  to  consider  the  purpose  of  this 
work,  namely  to.  state  tho  results  of  this  report  in  aerodynamic 
nomenclature.  Since  this  report  is  concerned  with  stability,  the  only 
contribution  of  the  aerodynamic  coefficients  is  how  they  appear  in  the 
basic  differential  equations.  This  means  that  in  order  to  obtain  the 
partner  of  K.  we  see  what  coefficient  appears  in  the  corresponding 
point  of  the"differential  equation  similar  to  (l)  which  is  based  on 
the  aerodynamic  force  system  (see  j,27j|  for  example).  By  this  tactic 
we  havet 

(C7)  ^ <V  4CmJ 

d q d 

to  when 

(08)  Kg->lA£§  (C„  ♦CH.). 

■ ' d'  q a 

• , - 

Notet  The  method  of  obtaining  (C7)  and  (C8)  is  not  too  desirable. 

It  would,  of  course,  be  more  satisfying  to  enlarge  the  ballistic  force 
system  so  that  there  would  exist  a one-to-one  correspondence.  (Efforts 
are  being  made  in  this  direction  at  the  present  time.)  It  also  she  Id 
be  noted  that  (C7)  follows  from  a comparison  of  the  homogeneous  equations. 
In  the  yaw  of  Vepose,  equation  (11),  should  be  replaced  by  - 1/h  c2S  CL, 


Since  the  major  function  of  Kg  is  its  contributions 
the  center*  of  mass  is  altered  we  havet 


By  use  of  (02),  (06),  (C?)  and  (08)  itf  is  now  possible  to  convert 
:•  our  symbols*  We  will  merely  tabulate  the  results.1  (K^  will  be  replaced 


Xjj  — Kjj  for  this  purpose.) 

H-^ij 

[\  * 2"°D  * V2 

(4>8  \ * S.-Q 

\*T<$  4 

M - It,*2 

*S  °«0 

•-/T? 

P»#  + °D  * *1 

M 

v) 

1* 


In  order  to  avoid  confusion  cl*  in  *5  (a*) 


will  be  replaced  by  y»* 
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